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Abstract 

We include new data in an updated analysis of helium in low metallicity extragalactic 
H II regions with the goal of deriving the primordial abundance of '^He (Yp). We show that 
the new observations of Izotov et al. (ITL) are consistent with previous data. However 
they should not be taken in isolation to determine Ip due to the lack of sufficiently low 
metallicity points. We use the extant data in a semi-empirical approach to bounding the 
size of possible systematic uncertainties in the determination of Y-p. Our best estimate 
for the primordial abundance of ^He assuming a linear relation between "^He and 0/H is 
Ip = 0.230 ± 0.003(stat) based on the subset of Hii regions with the lowest metallicity; 
for our full data set we find Ip = 0.234±0.002(stat). Both values are entirely consistent 
with our previous results. We discuss the implications of these values for standard big 
bang nucleosynthesis (SBBN), particularly in the context of recent measurements of 
deuterium in high redshift, low metallicity QSO absorption-line systems. 



1 Introduction 



Determining the primordial abundances of the Ught elements D, '^He, ^He, and ^Li is 
crucial for testing the standard model of big bang nucleosynthesis (SBBN) and in us- 
ing SBBN to set constraints on cosmology (e.g., the baryon density) and on particle 
properties (e.g., the number of light degrees of freedom contributing to the energy den- 
sity of the early universe (cf. Walker et al. 1991)). Each of the light nuclides poses 
different problems in the quest for its primordial abundance. For example, deuterium 
is destroyed during the pre-main sequence evolution of stars. As such, the relation be- 
tween its observed and primordial abundances depends on the amount of material which 
has been cycled through stars (i.e., on Galactic chemical evolution). Uncertainties in 
Galactic evolution currently dominate any attempt to infer the primordial abundance 
of deuterium from local ("here and now") observations although such data do provide a 
lower bound to its primordial value. 

The situation for ^He is much worse (i.e., even more subject to uncertainties in 
stellar and Galactic evolution) since stars of all masses destroy some - but not all - of 
their prestellar ^He and low mass stars are expected to be net producers of ^He. Thus, 
to infer the primordial abundance of ^He from observations "here and now" requires 
that the balance between destruction, survival and new production be known to an 
accuracy that eludes us at present. Since ''Li is observed in extremely metal-poor stars. 
Galactic chemical evolution plays virtually no role in using the data to infer its primordial 
abundance. However, since the very metal-poor stars are also very old, they have had 
a long time to modify their surface abundance of ^Li from its nearly primordial initial 
value. Uncertain corrections for stellar evolution dominate the uncertainty in the derived 
primordial abundance of lithium. 

Since "^He is the most abundant nuclide in the universe (after hydrogen) it may be 
observed throughout the universe (not only "here and now") more accurately than any 
of the other nuclides. For individual, careful observations, the abundance of ''He may 
be determined to 5% or better (0.012 in Y, the ^He mass fraction). Since its primordial 
abundance is expected to be quite large (~ 25% by mass) and ^He has been observed 
in the less processed, metal-poor, extragalactic Hii regions, the extrapolation to pri- 
mordial is minimal (differences in Y ^ 0.002 - 0.004). However, in the context of 
SBBN, the predicted primordial abundance of "'He is rather insensitive to (only loga- 



2 



rithmically dependent on) the one free SBBN parameter - the nucleon-to-photon ratio 
Tj (t) = riN/njirjio = lO^^rj). For this reason it is necessary to estabhsh the primor- 
dial abundance of ^He very accurately. Particular attention must be paid to possible 
systematic uncertainties as well as to the usual statistical errors. 

Although, compared to the other hght isotopes, the primordial ^He abundance is 
rather insensitive to chemical and stellar evolution, corrections for such effects are not 
entirely absent. In the course of their evolution stars burn hydrogen to helium and when 
they die they return this processed material (containing new ^He along with heavy el- 
ements ("metals") such as C, N, O) to the interstellar medium (ISM) polluting the 
primordial ^He. To minimize the contribution from stellar-produced ^He, the best ob- 
servational targets are those regions whose heavy element abundances are low, suggest- 
ing the least contamination from stellar and galactic chemical evolution. This has led 
virtually all investigators to concentrate on the low metallicity, extragalactic Hii re- 
gions (Searle & Sargent 1971; Peimbert & Torres- Peimbert 1974; Lequeux et al. 1979; 
Kunth & Sargent 1983; Torres-Peimbert et al. 1989; Pagel et al. 1992, hereafter PSTE; 
Skillman & Kennicutt 1993; Skillman et al. 1994, 1996, hereafter S; Izotov, Thuan, 
& Lipovetsky 1994, 1996 hereafter ITL). However, since in even the lowest metallic- 
ity regions observed (with metallicity as low as 1/50 of solar) some ^He was produced 
along with the heavy elements, Peimbert & Torres-Peimbert (1974) proposed to corre- 
late loBS with metallicity and to extrapolate to zero metallicity in order to infer Yp. 
Since the heavy element mass fraction, Z, is not observed, the abundances of oxygen 
and/or nitrogen have usually served as surrogates (e.g., Z 20(O/H)). 

A previous analysis (Ohve & Steigman 1995, hereafter OS), considered the ^He, 
0/H and N/H data of PTSE and of S, for 49 separate, low metaUicity, extragalactic 
Hii regions. In that analysis the correlation between 0/H and N/H was explored, and 
it was concluded that the nitrogen observed in these objects is dominated by a primary 
contribution (i.e., N/H increasing linearly with 0/H) with a small but not entirely 
insignificant secondary component. The fits of Y versus 0/H or Y versus N/H were not 
sensitive to the inclusion or exclusion of galaxies which show Wolf-Rayet features or of 
galaxies which deviated from the mean N/H versus 0/H relation. There was, however, 
a small difference between the values of Yp derived from the fits of Y versus 0/H and Y 
versus N/H. Given that there are both primary and secondary contributions to nitrogen, 
this is not unexpected (Fields 1996); we will return to this issue below. Overall, it was 
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found that the data were well described by a hnear fit of y to 0/H with an intercept at 
zero metallicity of Ip = 0.232 ± 0.003. In addition to the above statistical error, various 
contributions were described which might lead to an overall systematic uncertainty of 
order ± 0.005 (see also PSTE). 

In the past year or two the ^He abundance, key to the consistency of SBBN, has 
come under great scrutiny. Most recently, ITL presented new data which they claimed 
provides evidence for Yp in excess of 0.24. In this paper, we consider the ITL data and 
ask if they are consistent with those of PSTE and S. We conclude that they are and 
we propose an explanation for the apparent contradiction. We then use all extant data 
(PSTE, S and ITL) to derive the current best estimates for the primordial abundance 
of ^He. We also attempt to use the data in a semi-empirical approach to estimating the 
size of the possible systematic uncertainty in Ip. 

2 The Old Data used in OS 

OS used the data of PSTE and S for 49 separate extragalactic H ii regions. In minimizing 
the extrapolation to zero metallicity, the lowest metallicity Hii regions play a crucial 
role. As a result a "first cut" was made in OS eliminating those Hii regions (albeit only 8 
out of 49) with N/H > 1.0 x 10"^ and 0/H > 1.5 x 10"^. Note that all of the Hii regions 
retained are metal-poor compared to the Sun where (0/H) q = 8.5 x 10""^ and (N/H) q 
= 1.1 X 10"^. Nevertheless, the OS "first cut" data set spans one order of magnitude in 
oxygen abundance (15 ^ 10^{O/H) ^ 150) and a factor of ~ 25 in nitrogen abundance 
(4 ^ 10'^{N/H) ^ 100). OS also considered an even more metal-poor subset ("second 
cut"), retaining the 21 (out of 41) Hii regions with 0/H < 8 x lO'^ ([0/H] < -1). 
This more restricted 2nd cut set still has a modest dynamical range in its oxygen and 
nitrogen abundances: 15 < lO^iO/R) < 80 and 4 < 10^(N/H) < 40. 

As mentioned above, OS investigated the correlation between N /H and O /H for these 
H II regions. Although the variation of nitrogen with oxygen is of interest for the study 
of chemical evolution, it must be emphasized that the evolution of the very low mass 
host galaxies of these extragalactic H ii regions is likely dominated by local - in space 
and in time - processes. Different H ii regions may be "caught" at different evolutionary 
epochs (e.g., just before or just after a starburst, shortly before or immediately after a 
supernova explosion, etc.). Overall, OS found a strong correlation between N and O, 
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and that at low metallicity nitrogen is predominantly primary (varying linearly with 
oxygen) with a small, but not insignificant, secondary component (proportional to the 
square of the oxygen abundance). This can be seen from a power law fit to the data, 
N/H oc (0/H)", where OS found a. — 1.31 ± 0.07. This behavior is confirmed with our 
enlarged data set, now containing 62 distinct extragalactic Hii regions (labelled set B 
below), for which we find a — 1.21 ± 0.06. Alternatively, the predominantly primary 
nature of nitrogen can be seen by fitting the data with a hnear N/0 versus 0/H relation 
for which we find: 

N/O = (0.023 ± 0.002) + (76 ± 19){0/H) (1) 

Note that the "primary" component dominates for 0/H < 3.1 x 10"^, and that for our 
entire "first cut" range the "secondary" to "primary" ratio varies from 5 - 50 %. This is 
in agreement with Pagel & Kazlauskas (1992) who concluded that "primary" nitrogen 
dominates at low metallicity. 

Pagel, Terlevich & Melnick (1986) noted that Hii regions which showed Wolf-Rayet 
spectral features often had larger abundances of both helium and nitrogen compared to 
H II regions with the same oxygen abundance but lacking such features. OS searched for 
such an effect but found no statistically significant correlation, so OS did not exclude 
any H ii regions with Wolf-Rayet features and neither will we in our analysis here. This 
conclusion is supported by ITL and by Kobulnicky & Skillman (1996). 

In the previous Y versus 0/H analysis, OS found for the 1st cut (2nd cut) set 
Yp = 0.232 ± 0.003 (0.229 ± 0.005). Since then the data on extragalactic Hii regions 
has increased significantly. Izotov et al. (1994) presented observations of 10 Hii regions, 
four of which overlap those in the S set used in OS. These new data were incorporated 
in the analysis of Ohve & Scully (1996) who found that Yp derived from Y versus 0/H 
for the expanded 1st cut set (now containing 47 H ii regions) is slightly higher (although 
within 1(7 of OS): Yp = 0.234 ± 0.003. 

3 The New Data of ITL 

In their most recent work Izotov et al. (1996) have data from observations of 28 new 
H II regions. Of these one region is included in their 1994 set and four others are contained 
in the PTSE set, including a reobservation of IZwl8. Thus we now have data for 78 
distinct extragalactic H ii regions (several of which have been observed by two or more 
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independent groups). In the following analyses we impose the same low metallicity 1st 
cut restriction as in OS, ehminating the same 8 regions from the PTSE set. Thus the 
largest, low metaUicity set available for analysis consists of 70 regions (set A). However, 
according to several criteria ITL exclude 10 regions (including IZwl8) from their analysis; 
here we accept their judgment and with one notable exception (see below) eliminate 
the same regions from consideration. Since eight of the ITL excluded regions are not 
contained in the PSTE or S sets, our data set is now reduced to 62 distinct H ii regions 
(set B). For our second cut (in this case at 0/H < 8.5 x 10~^; [0/H] < -1) we have 32 
distinct Hii regions (set C). However, before proceeding it is necessary to consider the 
relation between the ITL data and those of PTSE and S. Are they consistent? 

One issue concerns the calculation of the statistical errors reported with the observed 
line ratios (corrected for reddening) . PSTE report errors derived from the total counts 
in the line and the continuum, and terms accounting for the sky subtraction and the read 
noise of the detector. Errors were not calculated for the reddening correction, the flat 
field correction, or the wavelength-dependence of the sensitivity (the "flux" cahbration), 
but care was taken such that these errors were of order, or smaller, than the calculated 
errors (see Simonson 1990). Skillman et al. (1994) include all of these terms in their 
uncertainties (see e.g., equation 2 in Skillman et al. 1994). 

ITL do not provide sufficient details to permit us to determine how the errors in 
their emission line ratios are calculated. But they do state that their spectra are "in 
excellent agreement" with those previously published in the literature. Frequently, their 
brighter lines are quoted with errors between 0.1 and 0.2% . This is to be contrasted 
with the analyses of PSTE and S where the minimum uncertainties usually lie in the 
range of 1 to 2% . It is interesting that ITL do comment that the residuals in their 
flux cahbration curve are " < 5% ." From this statement it appears possible that ITL 
might have underestimated the uncertainties in some of their reported line ratios by 
at least a factor of 25 (this may be a lower limit which could increase if other errors, 
e.g., the reddening correction, are accounted for). The possibihty that ITL may have 
significantly underestimated their observational errors is supported by, for example, the 
5(7 change in the reported oxygen abundance for 0940+544N from 7.37 ± 0.02 in ITL94 
to 7.48 ± 0.01 in ITL96. 

There are further reasons to suspect that ITL may have underestimated their errors. 
For example, it is not possible to know relative line ratios with an accuracy better than 
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the calibration of the telescope/spectrograph/detector combination. In assembhng a set 
of spectrophotometric standard stars for use with the HST, Oke (1990) used CCDs that 
were measured to be hnear to within 0.2% and found that standard star measurements 
were repeatable to "about 1% over most of the spectral range and a little larger in the 
ultraviolet and near infrared." Thus, it would be prudent to adopt 1% as a reasonable 
lower limit to the error on any measured emission line ratio. 

It is beyond the scope of this paper to devise our own corrections to the error esti- 
mates of ITL. Indeed, in some cases the hne ratio errors may not be very far off (e.g., 
for the weakest lines, the quoted errors become comparable to and even surpass the 
uncertainty in the flux calibration). Thus, for our analysis we simply adopt the ITL 
reported uncertainties. 

ITL have adopted an approach of using the data for four Hel lines (A4471, A5876, 
A6678, A7065) in a self-consistent analysis whose goal it is to determine simultaneously 
the recombination and the coUisional excitation contributions to the observed emissiv- 
ities. By insisting that the line ratios have their recombination values after correction 
for coUisions, ITL determine the electron densities self-consistently. The virtue of this 
approach is that it avoids the use of uncertain electron densities determined indirectly 
from [SII]. The problem with this approach is its rehance on the A7065 line which, al- 
though sensitive to coUisional excitation (albeit with an uncertain collision strength), is 
well-known (Robbins 1968; PSTE; G. Ferland, Private Communication) to be subject 
to fluorescence. The observed A7065 line strengths may well be providing a measure 
of the optical depth through the H ii regions rather than of the effect of coUisional ex- 
citation (although ITL argue to the contrary). Unquantified radiative transfer effects, 
complicated by unknown Hii region geometry, dust/gas, etc., may introduce large un- 
certainties in the ITL approach which call into question the efficacy of their reliance on 
this line. Further, their approach requires that ITL have good data for all four hues and 
this forces them to reject otherwise good observations of Hii regions when they have 
sufficiently accurate data for only two or three of the four hues. In contrast to the ITL 
approach, neither PSTE nor S use the A7065 line in their analyses and Peimbert (1996) 
notes that for the relatively low electron densities common to H ii regions the coUisonal 
correction is usually quite small. 

ITL use their method to analyze their data in a number of different ways. They have 
compared the He abundances derived from the observed emission line strengths using the 
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independent sets of recombination line emissivities by Brocklehurst (1972; hereafter B72) 
and Smits (1996; hereafter S96). For AA 5876, 4471, and 6678, the new S96 emissivities 
are in good agreement with those of B72 (see Table 3 of S96 for a comparison). Over 
the relevant range of electron temperature and density, the B72 and S96 emissivities 
for A5876 and A6678 agree to within one percent and for A4471 to within two percent. 
There is roughly a 40% difference for A7065 where B72 was in error (see discussion in 
Smits 1991a,b). Given the error in the B72 A7065 emissivities, it does not make sense 
for ITL to use the B72 emissivities in concert with their method. 

ITL also use two different sets of coUisional excitation rates to correct for the con- 
tribution of emission from coUisional excitation from the metastable 2^S level of He I. 
Clegg (1987; hereafter C) calculated these rates from the 19-state (up to n = 4) R-matrix 
computation by Berrington & Kingston (1987). Kingdon & Ferland (1995; hereafter KF) 
have calculated new rates, based on the 29-state (up to n= 5) computation of Sawey 
& Berrington (1993). Figure 1 compares the results of the C and KF calculations for 
the relative rates of coUisionally excited emission to recombination emission (C/R) for 
the four He lines used by ITL. Note the excellent agreement for A5876 and A6678. The 
small change in the rate for A4471 is due to both a change in the rate to the n = 4 level, 
and to the addition by KF of rates to two n = 5 levels. The large change in the A7065 
C/R value is mainly due to the difference between the recombination rates. C used B72 
for these rates and KF used S96 for these rates. Since the emissivities are implicit in 
the C/R calculations, it makes no sense to combine C with S96 nor KF with B72 as ITL 
have done. While the differences for A5876 and A6678 will be negligible, the differences 
for A4471 will be significant, and for A7065 very large. 

Therefore, when using the ITL analysis, the recombination emissivities must be 
restricted to those from S96 while any analysis which avoids the A7065 line may use either 
B72 or S96. Since previous analyses (PSTE; S) have avoided this hne (and have employed 
B72 /C) , it is interesting to compare the Y values derived from the ITL observations using 
S96/KF (their "best" combination) to those derived using B72 for the recombination 
emissivities, C for the collision strengths, along with the electron densities determined 
from [SII] . The ITL data reveal these differences to be quite small, in fact a weighted 
average of difference of the ^He mass fractions, ( Y(S96/KF) - Y(B72/C/SII) ) = -0.003 
which is much less than the typical errors in the individual Y determinations. Indeed, 
most of this difference is traceable to the fa 1% differences between the KF and C coUision 
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strengths which, however, lead to no significant difference (i.e., < 0.2%) when ITL 
evaluate Yp using S96/KF or S96/C (see Table 7 of ITL 1996). Within the uncertainties, 
Y(S96/KF) = Y(B72/SII). This agreement suggests that the analyses of PSTE and of 
S have not been biased by their rehance on B72, and that their results using B72/C/SII 
should be directly comparable to those of ITL. 

In contrast, ITL found from linear regressions of Y on 0/H using their own data 
Yp(S96/KF) = 0.243 ± 0.004 whereas OS, using PSTE and S data, found Yp(B72/SII) = 
0.232 ± 0.003. Given the relatively small statistical uncertainties such a large difference 
suggests an inconsistency between PSTE and S on the one hand and ITL on the other. 
Some have interpreted this apparent discrepancy as an indicator of the true size of the 
systematic errors in Yp determinations and have embraced the larger ITL values as a 
better probe of Yp. ITL apparently believe the inconsistency is real and they claim it 
is traceable to the use by PSTE and S of the older B72 emissivities. But we have just 
demonstrated, using the ITL data, that this cannot be the case; within the errors for all 
ITL Hii regions: Y(S96/KF) = Y(S96/C) = Y(B72/SII). Therefore, we must be able 
to compare the value of Yp derived from the ITL observations using the B72/C/SII 
combination with that found by OS using the PSTE and S data. For the data from the 
ITL preferred set of 27 Hii regions, a linear regression of Y on 0/H yields Yp = 0.241 
± 0.004, still quite high compared to the OS result. 

Why, then, do ITL find such a large value for Yp (0.241) compared to that found by 
OS from the data of PSTE and of S (0.232)? We believe it is due to the absence of the 
very lowest metallicity Hii regions from the ITL set after they have excluded selected 
regions from their analysis. In Figure 2 we show (with the error bars suppressed) the Y 
versus 0/H data used in OS (open circles) along with the new ITL data (filled circles). 
The crossed circles are the 10 Hii regions that ITL excluded from their analyses. Note 
that where there is overlap in 0/H, the ITL Y values are intermingled with those from 
PSTE and S. Indeed, there are six Hii regions in common between ITL and OS (not 
counting those ITL regions, shown in Figure 2, which ITL excluded from their fits) with 
four having higher and two with lower Y values than their PSTE or S counterparts. 
It is indeed surprising that ITL and OS find such significantly different values for Yp. 
However note that the ITL data set does not extend to as low an oxygen abundance 
as the set employed by OS. This more limited range in metaUicity for the ITL data set 
gives them less leverage in determining the slope of the Y versus 0/H relation. Indeed, 
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in the Y versus 0/H fit of the ITL data set, the slope is found to be 64 ± 48. This 
is reminiscent of the earher work of Kunth & Sargent (1983) whose hmited metaUicity 
range also led to an indeterminate slope. Kunth & Sargent (1983) made the appropriate 
choice based on their data and took a weighted mean of their data to determine Yp, 
finding a high value of 0.245 ± 0.003 which is very similar to those found by ITL. Indeed, 
when ITL used S96/KF or S96/C they found slopes consistent with zero which suggests 
that their Yp estimates are effectively weighted means. 

To test this hypothesis, we have refit the Y versus 0/H relation for the OS set 
excluding the four lowest 0/H points. For this modified "1st cut" OS set of 37 H ii regions 
we find Yp = 0.237 ±0.004, significantly higher than the previous result for the full "1st 
cut" set (0.232 ± 0.003). This refiects the high weights in the fits carried by the lowest 
metaUicity Hii regions which tend to have low Y values with small uncertainties. It is 
this modified OS value of 0.237 which should be compared to the ITL result of 0.241. 
Within the statistical uncertainty they are entirely consistent. To explore this further, 
we have employed the modified "1st cut" fit to describe the ITL data. The reduced 

for this fit is 0.50 (to be compared with 0.44 for ITL's own fit). Based on the F- 
test (Bevington & Robinson 1992) there is a 39% chance that their data is drawn from 
a distribution described by our fit. In addition, we have used the modified set of 37 
Hii regions in a "statistical bootstrap" of 40,000 runs (see Olive & Scully 1995) and 
we found that Yp exceeded 0.241 13% of the time. This is shown in Figure 3. These 
tests lead us to conclude that the new ITL data do not differ statistically from the older 
PSTE and S data used by OS. 

With this as justification, we proceed to analyze the combined data of PSTE, S 

and ITL. In this analysis we adopt the ITL data derived using the electron densities 

determined from SII so that we may have an internally consistent data setQ. Although we 

avoid using the ten Hii regions discarded by ITL, eight of them are rather insignificant 

in the sense that since their abundances have large uncertainties (and another is at 

intermediate metaUicity), they would have low weight in our fits and not much influence 

on our derived value of Ip. The one exception is IZwl8 which provides the lowest 

metaUicity point and for which ITL have good data. ITL exclude IZwl8 from their 

*We have used the values of Y and ay given by ITL. We note that the error in Y was not statistically 
propagated from that in the abundance by number, y. Thus the quoted errors are somewhat larger (by 
~ 30%) than they should be. We did not correct for this. 
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analysis because the A5578 line is subject to absorption by interstellar sodium. Not 
being able to use all four Hel lines in their self-consistent approach, ITL eliminate 
IZwl8 from their analysis. Izotov (Private Communication) has kindly provided us with 
the ^He abundance derived from their data using the A6678 hne and the [SII] electron 
density and we use this information to include this data point in our subsequent analysis. 
Note that when IZwlS is included with the other 27 ITL data points, their intercept 
drops to 0.239 ±0.004 (the effect is minor since the error is relatively large - in OS the 
error in IZwlS is diminished due to multiple as well as high quality observations). If we 
make our second cut to the ITL data including IZwlS, the intercept drops to 0.231 ± 
0.006. This result is now completely consistent with the OS value of 0.229 ±0.005. 

4 Results 

4.1 Yp From The Helium-Metallicity Correlation 

We adopt several approaches to using the H ii region data (set B with 62 separate H ii re- 
gions and the lower metallicity set C with 32 regions) to infer the primordial abundance 
of ^He, Yp. Since primordial hehum has been contaminated with the debris of stellar 
ejecta, the most common approach has been to use the metallicity information to probe 
the correlation of Y with Z (either 0/H or N/H) and to extrapolate this empirical rela- 
tion to zero metallicity to find Yp (Peimbert & Torres- Peimbert 1974). With increasing 
numbers of very low metallicity Hii regions, this extrapolation is quantitatively quite 
small (Ay 0.002 - 0.004). In OS we showed that the extant data do, indeed, justify 
a positive correlation between Y and 0/H (N/H). To explore this from a somewhat 
different perspective, consider the following: For set B we have computed the weighted 
mean of the hehum abundances, {Y) and in Figure 4 we plot the residuals, Y —(Y) 
as a function of the oxygen abundance. At low metallicity almost all the residuals are 
negative, while the positive residuals appear only at higher metallicity. Thus a one 
parameter fit to the Y versus 0/H data fails to account for the clear helium - oxygen 
correlation and is a poor fit to the data. Therefore we next try to fit the data with linear 
Y versus 0/H (N/H) regressions. These two-parameter fits describe the data very well 
(see Table 1). Note that if instead of our new "first cut" set B we had used the data 
for all 70 independent Hii regions (set A), there would be no difference in our derived 
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value of Yp. Similarly, there is no difference exceeding 0.001 in Ip if we exclude the ITL 
value (Izotov, Private Communication) for IZwl8 from our fits. 

Prom Table 1 we notice that the values of Yp derived from the Y versus N/H rela- 
tions are systematically higher (but only by ^ la) than those inferred from Y versus 
0/H. This effect is also present in OS and is entirely to be expected on the basis of the 
primary /secondary origin of nitrogen (Pields 1996). As Pields (1996) shows, the pri- 
mary/secondary origin for nitrogen, compared to the primary origin for oxygen implies 
that when Yp is derived from a linear correlation with N/H the result will exceed the 
"true" value derived from the linear Y versus 0/H relation. The quantitative difference 
between the two regressions will depend on the details of chemical evolution models as 
well as on the observed correlation of N with O and will be explored in future work 
(Pields, Olive & Steigman 1996). Por this reason we adopt for our estimates of Ip those 
values derived from the Y versus 0/H regressions for the B(C) sets (Table 1), 

Yp = 0.234 ± 0.002(0.230 ± 0.003) (2) 

for which Y^"" < 0.239(0.237). If we assume (PTSE) that the metaUicity Z and 0/H may 
be related by Z ^ 20{O/H), the set B(C) Y versus 0/H fits imply: AY/AZ 6(12), 
consistent with the steep slopes found by PSTE, Olive, Steigman, & Walker (1991) and 
OS. 



Table 1: Linear Fits for Y versus. O/H 



Set 


# Regions 


xVdof 


Yp 




10 slope 




0/H: B 


62 


0.58 


.234 ± 


.002 


1.25 ±0.27 


0.239 


0/H: C 


32 


0.59 


.230 ± 


.003 


2.31 ±0.65 


0.237 


N/H: B 


62 


0.61 


.237 ± 


.002 


25.7 ±5.9 


0.241 


N/H: C 


32 


0.63 


.232 ± 


.003 


61.5 ±17.5 


0.238 



4.2 Yp From A Few Good Hii Regions 

Por the fits described above, the extrapolation from the lowest metaUicity Hii regions 
to zero metaUicity is minimal (Ay 0.002 — 0.004). Nonetheless, it is true that any 
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extrapolation to zero metallicity could be avoided since the helium abundance inferred 
from the observations of any one Hii region (with non-zero metallicity) should provide 
an upper limit to Yp. For a very low metallicity Hii region such an upper limit may 
even provide a reasonable estimate of Yp. In this context, IZwl8, the most metal-poor 
Hii region, is an ideal candidate (Kunth et al. 1994) since it has been the subject of 
careful study by several independent groups (PSTE, Skillman & Kennicutt (1993) and 
ITL). A weighted mean of the five observations of the two separate knots in IZwl8 yields, 

Y(IZwl8) = 0.230 ± 0.004 (3) 

with a 2a upper bound of 0.237. In terms of statistical accuracy this result is fully 
competitive with the value of Yp derived in the previous section from 62(32) Hii regions. 

Of course, it should be kept in mind that the abundance inferred for any one H ii re- 
gion might be anomalous. Therefore the value of Y derived from the average of several 
Hii regions is also of interest. In such an analysis, as more regions are included, the 
mean value (weighted) of Y will increase, but if the errors are statistical, the error in 
the mean will decrease. As a result, for Hii regions the one-(or two-)cr upper bound 
to Y will first decrease with N, then level off and, as A" is further increased, it will 
eventually increase monotonically. This behavior is seen in Figure 5 where we show the 
weighted means, and the 2a bounds to the weighted means of Y derived from the A" 
lowest helium abundance Hii regions. Note that for 2 < A" < 13, the mean varies from 
0.229 to 0.231 while for 2 < A^ < 14, (Y) < 0.236 (2ct). It is not unreasonable to infer 
from these results that, 

Yp < 0.230 ± 0.003 (4) 

with, Yp^ < 0.236. If, instead, wc take the weighted means of the regions with the lowest 
values of 0/H, wc obtain a similar result. This illustrates the potentially great value of 
very careful analyses of a handful of the lowest metallicity (lowest Y) H ii regions. 

4.3 The Systematic Uncertainty In Ip: A Semi-Empirical Ap- 
proach 

Many observers have identified numerous sources of uncertainty affecting HII region 
helium abundance determinations (see, e.g., Davidson & Kinman 1985; Dinerstein & 
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Shields 1986; PSTE; S; ITL; Peimbert 1996). Peimbert (1996) divides the errors asso- 
ciated with the determination of the primordial helium abundance into three groups: 
(I) errors in the determination of the line ratios; (II) errors in the interpretation of the 
line ratios; and (III) errors in the extrapolation to Z — 0. Here we try to infer a rea- 
sonable estimate of the systematic uncertainty in Yp by inspecting the various possible 
systematic effects in each group. 

In Group I the errors in determining line ratios can be attributed to measurements 
of the line ratios (including signal-to-noise in the line and sky subtraction), detector 
calibration, reddening corrections, and lack of corrections for possible underlying stellar 
absorption. All of the effects in Group I have been discussed in detail in Skillman et al. 
(1994) and in previous studies. To summarize, if detectors which have been tested for 
linearity (CCDs) are used, if several standard stars which have been previously observed 
with hnear detectors (preferably the HST standards of Oke) are observed, if the targets 
are restricted to those objects of high excitation and high Balmer line equivalent width, 
and if one accumulates in excess of 10,000 photons in each of the He lines used, then it 
is possible to achieve an accuracy of 2% in the relevant He/H hne ratios. Then, of all 
the effects described above, only unaccounted-for underlying stellar absorption would 
cause a systematic error, leading to an underestimate of the He abundance. However, 
the presence or absence of this effect can be probed by measuring different He lines of 
different equivalent widths. The general agreement between the different lines, in those 
cases with careful tests, indicates that the effect of underlying stellar absorption is of 
order 1% or less. 

In Group II the errors in the interpretation of line ratios can be attributed to correct- 
ing for the presence of neutral He, variations in temperature structure ("temperature 
fluctuations" ) , the accuracy of the atomic data, the correction for the coUisional excita- 
tion of Hel lines (primarily from the meta-stable triplet 2S level) , correction for radiative 
transfer effects, and correction for coUisional excitation of the HI hues (from the ground 
state). Taken in order: 

(i) In principle, the presence of neutral He would systematically lower the observed 
He abundance. However, none of the tests performed so far have found any evidence of 
neutral He (see Skillman et al. 1994), and photoionization models indicate that it is not 
hkely to be a problem for the objects included in these studies (see also the discussion 
in Vilchez & Pagel (1988) and in PSTE). 
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(ii) If there are variations in the electron temperature in the gas, the heavy element 
abundances derived from the coUisionally excited lines would be underestimated and the 
He abundances would be shghtly overestimated (due to the weak dependence of the He 
hues on the electron temperature, the effect is 1.5% for A5876 and A6678 for an error of 
lOOOK at 15,000K and about half of that for A4471). Temperature variations appear to 
be much more likely at higher metallicities, but if supernovae are an important heating 
source in the Hii regions, then temperature fluctuations may be important (Skillman 
1995, Peimbert et al. 1991). 

(iii) The comparison of Smits (1996) with Brocklehurst (1972) would indicate that the 
atomic data and calculations of the recombination emissivities for the He lines of interest 
are good to better than 1%. Note, however, that there are much larger differences in 
the infrared transitions, and more work (both theoretical and observational) is desirable 
in this area. 

(iv) The recent work by Kingdon & Ferland (1995) gives us confidence that we are 
able to correct accurately for coUisional excitation of Hel. These corrections are usually 
of order 1-3%. Not correcting would systematically overestimate the He abundance. 
The main problem here is to determine the density sufficiently accurately. ITL have 
argued that densities derived from [S II] emission lines are not appropriate. Since the 
different He I emission lines have different density dependences, by using several lines 
it is possible to solve for the density, and thus the correction (this is essentially how 
ITL propose to solve for the electron density). In general, the densities derived in this 
manner by ITL agree well with those derived using the [S II] hues. 

(v) Since the work of Robbins (1968), it has been generally agreed that radiative 
transfer effects are unimportantly small for the bright lines which are used to derive the 
He abundances, particularly the singlet lines. Based on the results of photoionization 
modeling, in which the effects of coUisional coupling of the singlets and triplets, radiative 
transfer effects, and coUisional excitation were all treated simultaneously, Sasselov & 
Goldwirth (1995) claimed that the He/H hne ratios lead to systematic underestimates. 
However, no He/H line ratios were presented in their paper. Until such effects on the 
He/H line ratios are identified and quantified, it seems reasonable to ignore this claim. 

(vi) Davidson & Kinman (1985) showed that at the high temperatures found in the 
lowest metallicity HII regions, coUisional excitation of the HI lines may be important. 
Skillman & Kennicutt (1993) showed that this effect, which is dependent on the neutral 
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hydrogen fraction, is not significant for neutral H fractions less than 0.0001. Straight- 
forward calculation of the photoionization balance in an Hii region usually results in 
neutral H fractions less than this. Photoionization codes often produce higher neutral 
H fractions, but this may be due to the approximations made in the treatment of the 
ionizing radiation field. This could be a very difficult uncertainty to pin down, since the 
geometry of the gas distribution relative to the ionizing source has a strong influence on 
the neutral H fraction. If this effect were important, it would result in an underestimate 
of the He abundance. Perhaps 2% represents a reasonable upper limit on the uncertainty 
of this effect. 

Group III errors are concerned with the extrapolation of the observed helium abun- 
dances to zero metaUicity. The "classical" approach (Peimbert & Torres-Peimbert 1974) 
has been to fit the observations with a linear Y versus Z relation (where the oxygen abun- 
dance usually serves as a surrogate for Z) and to extrapolate to Z = 0. For observations 
of low metaUicity H ii regions this linear fit may be thought of as the lowest order con- 
tribution to a more general Y(Z) relation. For the set B(C), this extrapolation from the 
lowest metaUicity data is quite smaU: AY = 0.002±0.001(0.003± 0.001). Since He/H 
is only expected to increase with Z, it is unlikely that this approach can systematically 
underestimate Yp. While this linear fit may yield an upper bound on Yp, it does not 
necessarily provide a lower bound. Indeed, as our second cut set (set C) shows, the Y 
versus Z slope appears to steepen at the very lowest metallicities. If the assumption of 
linearity is relaxed, then, in principle, Yp can be significantly lower. 

Given the different production sites of He and O, He/H may not be expected to 
track 0/H well. Indeed, Steigman, Gallagher & Schramm (1989) suggested that hehum 
may correlate better with nitrogen and/or carbon than with oxygen. The observation 
by Pagel et al. (1992) that the dispersion in the He vs. N regression is less than that of 
the He vs. O regression lends some support to this expectation. However, the observed 
hnearity of Y with Z (where Z fa 20(0/ H)) over more than a decade in Z (e.g., OS) 
may reflect a balance between losses due to galactic winds (most important for low 
mass, low metaUicity systems) and the metaUicity dependence of the yields (O yield 
decreasing with increasing metaUicity due to the increasing importance of stellar winds 
in the massive stars). Future work on accurate relative abundances to search for a 
"second parameter" in the Y vs. Z relationship would be of great value. 

Since the data are entirely consistent with a linear Y vs. Z relation, the uncertainty 
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in the intercept should be of order the uncertainties in the best measured points at 
low metallicities (i.e., 2 - 3%). Calculating the uncertainty in the intercept depends on 
knowledge of the source of the scatter in Y at a given Z, which presently is dominated 
by measurement uncertainties. 

In the analysis presented here we have, to some extent, avoided the issue of the 
extrapolation to zero metallicity by considering the helium abundance determined from 
the best observed Hii regions (Yp < Yobs)- From five independent observations of 
IZwl8 we found, Yp < 0.230±0.004; from the 13-14 Hii regions with the lowest hehum 
abundances we derived, Yp < 0.230±0.003. 

Finally, there remains the important question of how to combine different sources 
of systematic uncertainties. Since the possible errors are not correlated, it makes no 
sense to add (linearly) all imaginable systematic errors to obtain an estimate of the 
total systematic error. Many potential error sources can be classified as unlikely, with 
bounds constrained observationally. Therefore it is even more unlikely that a single data 
point (let alone all of them) would suffer from more than one of the potential systematic 
errors at the amplitude of the observationally constrained limits. The salient point is 
that all imaginable systematic errors appear to be limited to about 2 percent or less. 
Thus, it seems reasonable to adopt as an estimate for the overall systematic error, 0.005, 
as proposed by PSTE. 

We can also attempt to exploit the data itself to provide a bound on the size of 
possible systematic sources of error. Although many of the sources of potential errors 
listed above might shift Y and/or 0/H in a systematic fashion, modifying the intercept 
and/or slope inferred from the Y versus 0/H relation, their variation from source to 
source and observer to observer (telescope/detector to telescope/detector) would also 
contribute to the overall dispersion of the data around the "true" Y versus 0/H relation. 
We have therefore taken our best fit for set B (see Table 1) and examined the residuals, 
Y — Y{B), as a function of 0/H. For the variance of the residuals we find 0.007; this 
should be compared to the error estimates for individual Hii region Y determinations 
which are, on average, 0.010. It appears that the observers have been generous, perhaps 
overly generous, in their error estimates as first pointed out by PTSE. Indeed, this was 
already suggested by the small values of the reduced seen in Table 1. We stress that 
since the variance of the residuals is small, there is no real statistical significance to the 
ITL claims that S96/KF improves the scatter. 
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As a further check on the stabihty of our Yp estimates and to constrain many of 
the sources of possible systematic error, we have performed 40,000 runs of a statistical 
bootstrap (Olive & Scully 1995) using all the data (set A) with, and without the error 
estimates (we thank Tim Beers for suggesting this test to us and Sean Scully for doing the 
runs) . In Figure 6 we show the resulting distributions for Yp (when errors are included) . 
In both cases the distributions are closely Gaussian with Yp = 0.234 ± 0.002 and 95% 
CL upper bounds ^ 0.238. This test suggests that unless ALL values of Y should 
be systematically shifted (e.g., due to inaccurate atomic data), 0.238 might provide a 
good upper bound to Yp, including systematic errors. However, to err on the side of 
conservatism, instead of an upper bound of 0.238 (including systematic uncertainties), 
we will adopt the set C(B) results in the subsequent discussion, 

Yp = 0.230 ± 0.003 ± 0.005(0.234 ± 0.002 ± 0.005) (5) 

where 0.003(0.002) represents the (Gaussian) statistical error (Fj'^ < 0.237(0.239)) and 
0.005 is a possible systematic offset in Yp (leading to "firm" 2a upper bounds to Yp of 
0.242 or 0.244). We will also explore the consequences of a larger value for AFsys. 

5 Discussion 

First let us ignore any possible systematic uncertainty in our adopted value of Yp (eq. 5) 
to identify the range in the nucleon abundance (r/io) which follows from SEEN (including 
uncertainties in the neutron lifetime and any relevant nuclear reaction rates; Hata et al. 
1995). For Yp = 0.234 ± 0.002, 

7710 = 1.8 ±0.3 (6) 

For an upper bound of Yp < 0.239(0.244), the corresponding 95% CL upper bound on 
rjiQ is 2.4(3.8). Furthermore, the very low value of ?7 in (|^), which is derived directly 
from ^He, corresponds to a very low universal density of baryons, 

fig/i^ ^ 0.007 ±0.001 (7) 

Even for Yp < 0.239(0.244), Qph"^ < 0.009(0.014). For the lower value of Yp = 0.230 ± 
0.003, 

7/10 = 1.4 ±0.3 (8) 
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In this case for Yp < 0.237(0.242), the upper bound on 7710 is 2.1(3.2). Consequently, 
the density of baryons is even lower, 

Ob/i^ = 0.005 ± 0.001 (9) 

and for Yp < 0.237(0.242), QbH'^ < 0.008(0.012). 

Although these estimates for 77 are based solely on ^He, concordance of SBBN with 
the observational data requires that the same value of rj yield acceptable abundances 
for the other light elements. For example, in a likelihood analysis based on ^He and ^Li 
(the latter inferred from the metal-poor halo stars). Fields & Olive (1996) and Fields et 
al. (1996) found consistency for a similarly low value of 77. Using Yp — 0.234 ± 0.003, 
they found r/io = 1.8 as the most likely value and a 95% CL range from 1.4 to 4.3. 
Repeating this ^He/^Li analysis for the lower value of Yp — 0.230 ± 0.003 we find for 
the best fit value 7710 = 1.7 and a 95% CL range from 1.3 to 4.0. 

For such low nucleon abundances there is consistency (Dar 1995) between the SBBN 
predictions and the primordial abundances not only of ^He and ^Li, but also with the 
deuterium abundance as determined from observations of some QSO absorption line 
systems (Songaila et al. 1994; Carswell et al. 1994; Rugers & Hogan 1996). The very 
high D abundance does pose a challenge to Galactic evolution models since it requires 
that the Galaxy has destroyed more than 90% of its initial deuterium. To do so while 
avoiding the overproduction of heavy elements may require the presence of supernovae 
driven Galactic winds (Scully eA, al. 1996). When the likelihood analysis is extended to 
include the high primordial D/H, Fields et al. (1996) find that the peak value of 7710 is 
1.7 with a 95% CL range between 1.5 and 2.4; we find that this range remains essentially 
the same for either Yp = 0.234 or 0.230. The corresponding value for the combination 
flBh'^ is now 0.006 with a 95% CL range between 0.005 and 0.009. 

Alternatively we may use either of the two recent determinations of primordial deu- 
terium in high redshift, low metallicity QSO absorption systems to pin down 77 and 
the corresponding range for Yp predicted by SBBN which we may then compare to our 
adopted value (range) for Yp. 

5.1 High-D 

If the high abundance of deuterium derived from the observations of Songaila et al. 
(1994), Carswell et al. (1994) and Rugers and Hogan (1996) is representative of the 
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primordial D abundance, then 1.3 < 7710 < 2.7 and 0.231 < Fsbbn < 0.239 (95% CL) 
(Hata et al. 1996). This range, as aheady seen above, is in excellent agreement with our 
adopted range for Yp (eq. 5) derived from the data and may be used to infer a restrictive 
upper bound to the effective number of equivalent light neutrinos (AA^,^ = A^,^ — 3). For 
a systematic offset to Ip of AYgys — 0(0.005), 

AN^ < 0.5(0.8) (10) 

Notice that if Yp > 0.239, > 3.0 would be required. The 95% CL upper hmit on the 
number of light degrees of freedom from the likehhood analysis of Fields et al. (1996) 
based on ^He and ^Li is (Olive & Thomas 1996) AA^^ < 1.0 for Yp = 0.234 ± 0.002 and 
AA^^ < 0.7 for Yp = 0.230 ± 0.003 (in both cases a^ys = 0.005 was assumed). 

5.2 Low-D 

In contrast, if the deuterium abundances derived for two different lines-of-sight from 
the data of Tytler, Fan & Buries (1996) and of Buries & lytler (1996) provide good 
estimates of the true primordial value, then 5.1 < r]io < 8.2 and 0.246 < Isbbn < 0.252 
(95% CL) (Hata et al. 1996). Consistency with SBBN (A^^ > 3.0) can only be recovered 
if systematic effects in deriving Yp from the data have led to an underestimate by an 
amount Al^yg > 0.009 (i.e., consistency would require that Yp > 0.246 compared to our 
upper bound of 0.237). Alternatively, if Ay^ys < 0.005, A^,. < 2.7. 

6 Summary 

For nearly two decades low metallicity, extragalactic Hii regions have been studied as 
a probe of the primordial abundance of helium. From observations of ten such re- 
gions Lequeux et al. (1979) derived Yp = 0.233 ± 0.005. Using four carefully studied 
Hii regions, Torres-Peimbert et al. (1989) found Yp = 0.230 ± 0.006, estabhshing the 
competitiveness of quality with quantity. On the basis of nineteen extragalactic H 11 re- 
gions PSTE inferred Yp = 0.228± 0.005, and building on this data set OS added the data 
of S to find for 41 (21) low metallicity Hii regions Yp = 0.232±0.003 (0.229 ± 0.005). 
In this paper we have considered the new data from ITL which, at first glance, seems to 
yield a much larger value for Yp. In contrast, we have found that the ITL data are fully 
consistent with those of PSTE and S and therefore we have combined these sets in an 
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analysis of 62 (32) Hii regions. Prom the Y versus 0/H correlation we derive (see Table 
1) Yp = 0.234 ± 0.002 (0.230 ± 0.003). We have also considered the five independent 
observations (PSTE; Skillman & Kennicutt 1993; ITL) of the two knots in IZwl8 to 
derive an upper bound to Ip of 0.230 ± 0.004, and in an extension of this approach to 
Yp we have considered the weighted mean (and the weighted mean plus 2-sigma) for the 
Hii regions with the lowest values of Y to derive Ip < 0.230 ± 0.003. These results have 
led us to adopt a "95% CL" upper bound to Yp of 0.237 (for set B with 62 Hii regions, 
this upper bound is 0.239). 

The availability of large numbers of carefully observed, low metallicity H ii regions 
has permitted estimates of Yp whose statistical uncertainties are very small 1%). 
However, there remains the possibility that in the process of using the observational 
data to derive the abundances, contamination by unacknowledged systematic errors has 
biased the inferred value of Yp. The observers have identified many potential sources of 
such systematic errors (Davidson & Kinman 1985; PSTE; S; ITL; Peimbert 1996) and, 
where possible, have designed their observing programs to minimize such uncertainties 
and/or to account for them. Here we have noted that many of the identified sources of 
potential systematic errors would vary from H ii region to H ii region and from observer 
(telescope/detector) to observer, introducing not only a systematic offset in the derived 
value of Yp, but also an accompanying dispersion in the Y versus 0/H relation. The 
very small values of the reduced for our fits suggest that the observers' error estimates 
may already account for some sources of systematic error. We have performed several 
tests confirming this and conclude that our determinations of Ip are robust in the 
absence of some yet to be identified systematic offset which shifts all the data uniformly. 
Nonetheless, in discussing the consequences of our derived value of Yp for cosmology and 
for particle physics we have allowed for a possible systematic offset AFgys = 0.005. 

For SBBN the low value we derive for Ip, consistent with previous results (Lequeux 
et al. 1979; Torres- Peimbert et al. 1989; PSTE; OS), imphes a low nucleon abundance 
but is entirely consistent with the inferred primordial abundances of ^Li and D (from the 
QSO absorbers studied by Rugers & Hogan 1996). Provided that the systematic error 
in Yp is not large, there is a meaningful constraint on the effective number of equivalent 
light neutrinos (Steigman, Schramm & Gunn 1977). If, instead, the low deuterium 
abundance inferred from the data of Tytler, Fan & Buries (1996) and of Buries & Tytler 
(1996) is the "true" primordial value, there is a challenge to SBBN unless AYgys is large 
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(so that Yp > 0.246). 



Acknowledgments 

We owe a debt of gratitude to R.C. Kennicutt, E. Terlevich and R.J. Terlevich for 
sharing their unpubhshed data with us. We thank S. ScuUy for running the statistical 
bootstrap and for preparing Figures 3 and 6. We also are pleased to thank T. Beers, D. 
Garnett, N. Hata, Y. Izotov, D. Kunth, B. Pagel, D. Schaerer, S. Viegas and T. Walker 
for very informative and valuable discussions. The work of KAO is supported in part by 
DOE grant DE-FG02-94ER-40823 and that of GS by DOE grant DE-AC02-76ER-01545. 
EDS is grateful for partial support from a NASA LTSARP grant No. NAGW-3189 and 
a Bush Sabbatical Fellowship from the Graduate School of the University of Minnesota 
and would especially like to thank the Max-Planck-Institute for Astrophysics for their 
hospitality during his sabbatical year. Parts of this work were carried out while KAO 
and GS were participants in the INT Workshop on Nucleosynthesis in Stars, Supernovae 
and the Universe and also while GS was a visitor at the Instituto Astronomico e Geofisico 
of the University of Sao Paulo (Brasil) and they are grateful for the hospitality provided. 



22 



References 



Berrington, J.A. & Kingston, A.E. 1987, J. Phys. B, 20, 6631 

Bevington, P.R. & Robinson, D.K. 1992 Data Reduction and Error Analysis for the 
Physical Sciences (New York: McGraw Hill) 

Brocklehurst, M. 1972, MNRAS, 157, 211 (B72) 

Buries, S. & Tytler, D. 1996, ApJ, 460, 584 

Carswell, R.F., Ranch, M., Weymann, R.J., Cooke, A.J. & Webb, J.K. 1994, MNRAS, 
268, LI 

Clegg, R.E.S. 1987, MNRAS, 229, 31P 
Dar, A. 1995, ApJ, 449, 550 

Davidson, K. & Kinman, T.D. 1985, ApJS, 58, 321 

Dinerstein, H.L. & Shields, G.A. 1986, ApJ, 311, 45 

Ferland, G.J. 1986, ApJ, 310, L67 

Fields, B.D. 1996, ApJ, 456, 478 

Fields, B.D. & Olive, K.A. 1996, Phys Lett B368, 103 

Fields, B.D., Kainulainen, K., Olive, K.A., & Thomas, D. 1996 New Astronomy, 1, 77 
Rata, N., Scherrer, R.J., Steigman, G., Thomas, D., & Walker, T.P. 1995, Phys. Rev. 
Lett., 75, 3977 

Hata, N., Steigman, G., Bludman, S. & Langacker, P. 1996, Phys. Rev. D, in press, 
astro-ph/9603087| 

Izotov, Y.I., Thuan, T.X., & Lipovetsky, V.A. 1994 ApJ 435, 647. (ITL) 
Izotov, Y.I., Thuan, T.X., & Lipovetsky, V.A. 1996, (ITL) 
Kingdon, J. & Ferland, G.J. 1995, ApJ, 442, 714 (KF) 
Kobulnicky, H.A. & Skillman, E.D. 1996, ApJ, 471, 211 

Kunth, D., Lequeux, J., Sargent, W.L.W., & Viallefond, F. 1994, A&A, 282, 709 

Kunth, D. & Sargent, W. 1983, ApJ, 273, 81 

Lequeux, J., et. al. 1979, A&A, 80, 155 

Oke, J. B. 1990, AJ, 99, 1621 

Olive, K.A. & Scully, S.T. 1996, IJMPA, 11, 409 



23 



Olive, K.A., Steigman, G. & Walker, T.P. 1991, ApJ, 380, LI 
Olive, K.A., & Steigman, G. 1995, ApJS, 97, 49 
Olive, K.A. & Thomas, D. 1996, preprint in preparation 
Pagel, B.E.J. & Kazlauskis, A. 1992, MNRAS, 256, 49P 

Pagel, B E.J., Simonson, E.A., Terlevicli, R.J. & Edmunds, M. 1992, MNRAS, 255, 325 
(PSTE) 

Pagel, B.E.J., Terlevicli, R.J. & Melnick, J. 1986, PASP, 98, 1005 (PTM) 

Peimbert, M. 1996, Rev. Mex. Astr. Astrofis., in press 

Peimbert, M. Sarmiento, A. & Fierro, J. 1991, PASP, 103, 815 

Peimbert, M. & Torres-Peimbert, S. 1974 Ap J, 193, 327 

Robbins, R.R. 1968, ApJ, 151, 511 

Rugers, M. & Hogan, C. 1996, ApJ, 259, LI 

Sasselov, D. & Goldwirth, D.S. 1995, ApJ, 444, L5 

Sawey, P.M.J. & Berrington, K.A. 1993, Atomic Data Nucl. Data Tables, 55, 81 



Scully, S.T., Casse, M., Olive, K.A., & Vangioni-Flam, E. 1996, ApJ, in press, astro- 



ph/9607106 



Searle, L. & Sargent, W.L.W. 1971, ApJ, 173, 25 
Simonson, E. A. 1990, Ph. D. thesis. University of Sussex 
Skillman, E. 1985, ApJ, 290, 449 
Skillman, E., & Kennicutt 1993, ApJ, 411, 655 (S) 

Skillman, E., Terlevich, R.J., Kennicutt, R.C., Garnett, D.R., & Terlevich, E. 1994, 

ApJ, 431, 172 (S) 
Skillman, E., et al. 1996, ApJ Lett (in preparation) (S) 
Smits, D.P. 1991a, MNRAS, 248, 193 
Smits, D.P. 1991b, MNRAS, 251, 316 
Smits, D.P. 1996, MNRAS, 278, 683 (S96) 

Songaila, A., Cowie, L.L., Hogan, C. & Rugers, M., 1994, Nature, 368, 599 

Steigman, G., Gallagher, J.S. & Schramm, D.N. 1989, Comments on Astrophysics, 
Vol. XIV, p. 97 



24 



Steigman, G., Schramm, D.N. & Gunn, J. 1977, Phys. Lett., B66, 202 
Torres-Peimbert, S., Peimbert, M., & Fierro, J. 1989, ApJ, 345, 186 
lytler, D., Fan, X.-M., and Buries, S. 1996, Nature, 381, 207 
Vilchez, J.M. & Pagel, B.E.J. 1988, MNRAS, 231, 257 

Walker, T. P., Steigman, G., Schramm, D. N., Olive, K. A., & Kang, H. 1991, ApJ, 376, 
51 (WSSOK) 



25 



Figure Captions 



Figure 1: A comparison of the calculations of the coUisionally excited emission 
to recombination emission rates (C/R) for four Hel recombination lines 
using the formulae of Clegg (1987) and Kingdon & Ferland (1995). The 
calculations were carried out for an electron density of 100, which is 
appropriate for most low metallicity HIT regions. The large difference seen 
for the A 7065 line arises because Clegg (1987) used the recombination 
coefficients from Brocklehurst (1972) which, for A 7065, have been shown 
to be in error by Smits (1991a,b). 

Figure 2: The hehum (Y) and oxygen (0/H) abundances for the extragalactic 
Hii regions of the 1st cut data sets (OS) from PTSE and S (open cir- 
cles), and from ITL (filled circles). Regions excluded by ITL are shown 
as crossed circles. Lines connect the same regions observed by different 
groups. 

Figure 3: The distribution of the error weighted determinations of Yp (open circles) 
from a 40,000 run statistical bootstrap using the modified OS data set 
(excluding the four points with lowest 0/H). The solid curve is the best 
fit gaussian. 

Figure 4: The residuals, Y — (Y), of the data (Y) compared to the weighted mean 
{(Y)) versus the oxygen abundance for the first cut set B. 

Figure 5: The running average (weighted means) of the "^He abundance, Y, for the 
first N (lowest Y) Hii regions. Also shown are the 2a bounds to the 
weighted means. 

Figure 6: The distribution of the error weighted determinations of Yp (open circles) 
from a 40,000 run statistical bootstrap using the full data set (A). The 
solid curve is the best fit gaussian. 
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